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Abstract 

Jiirgen Ehlers developed frame theory to better understand the relationship between 
general relativity and Newtonian gravity. Frame theory contains a parameter A, which 
can be thought of as 1/c 2 , where c is the speed of light. By construction, frame theory 
is equivalent to general relativity for A > 0, and reduces to Newtonian gravity for 
A = 0. Moreover, by setting e = vA, frame theory provides a framework to study 
the Newtonian limit e \ (i.e. c — ► oo). A number of ideas relating to frame theory 
that were introduced by Jiirgen have subsequently found important applications to 
the rigorous study of both the Newtonian limit and post-Newtonian expansions. In 
this article, we review frame theory and discuss, in a non-technical fashion, some of 
the rigorous results on the Newtonian limit and post-Newtonian expansions that have 
followed from Jiirgen's work. 



1 Introduction 

Throughout his scientific career, Jiirgen was interested in conceptual questions of physics, 
and in particular, the relationship between different theories that describe the same phe- 
nomena in nature. A prime example of this is the solar system. It can be, to high accuracy, 

*Todd. Oliynyk@sci.monash.edu.au 
^ B cmd . S chmidt @ aei . mpg.de 



1 



described using either Newton's or Einstein's theory of gravity. The question then becomes 
how can we understand this relationship between the two theories in spite of their quite 
distinct formulations. 

To begin to answer this question, Jiirgen developed, extending ideas of Cartan [3j and 
Friedrichs [TJ, a theory, called frame theory, which contains both Einstein's and Newton's 
theories of gravity. Frame theory contains a parameter A, which can be thought of as 1/c 2 
where c is the speed of light. By construction, frame theory coincides with Einstein's theory 
for A > 0, and reduces to Newton's theory when A = 0. In this way, frame theory provides 
a way of understanding the relationship between the laws, or, in other words, the equations 
of the two theories. 

In addition to developing frame theory, Jiirgen made another important contribution 
which he never published. He discovered (see Figure [TJ at the end of the article) a way of 
writing Einstein's field equations containing the parameter e = y/\ for which the partial 
differential equations of Einstein's theory (e > 0) become the equations of Newton's theory 
in the limit e \ 0. Following standard terminology, we will refer to the limit e \ as the 
Newtonian limit. 

Although, Jiirgen's formulation provided a satisfying method for addressing the relation- 
ship between the two theories at the level of equations, the problem of understanding the 
relationship between the solutions remained open. The relationship between the solutions is 
what one really wants to know as it is the solutions that describe the physical phenomena, 
and not the equations themselves. A first result in this direction was obtained by Lotter- 
moser [16]. Using Jiirgen's formulation, Lottermoser showed that there exists one-parameter 
families of solutions depending on e to the gravitational constraint equations that converge 
in the limit e \ to Newtonian initial data. This result was a necessary first step in the 
analysis of the limit e \ since it established the existence of reasonable, relativistic initial 
data that depends on e in the proper fashion. However, Lottermoser did not answer the 
existence question for relativistic solutions with a Newtonian limit. The first general result 
showing the existence of a Newtonian limit for any matter model was by Rendall for Vlasov 
matter. In [26], Rendall rigorously proved that Einstein- Vlasov system has a large class of 
solutions that have a well defined Newtonian limit. 

Another important existence result for the Newtonian limit, this time in the stationary 
setting, was obtained by Heilig [9J- Heilig proved the existence of stationary, rotating, axi- 
symmetric solutions of the Einstein-Euler equations using Jiirgen's formulation to perturb 
away from a Newtonian solution at e = to obtain a fully relativistic solution for e > 0. 
This paper still represents the only known existence proof for stationary, fully relativistic, 
perfect fluid solutions. Heilig also established that these solutions depend analytically on e, 
and therefore, admit convergent post-Newtonian expansions. Later, Heilig's methods were 
adapted by one of the authors to prove existence of static solutions for the Einstein field 
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equations coupled to Yang-Mills fields [TSJ, QIJ] . 

In [20], one of the authors adapted the variables introduced by Jiirgen to prove the 
existence of a wide class of solutions to the Einstein-Euler equations that have a Newtonian 
limit. We note that this formulation was different from Rendall's both in terms of the 
variable used to represent the gravitational field, and the choice of gauge. The methods used 
in [20] were subsequently generalized to prove the existence of a large class of solutions to the 
Einstein-Euler equations that admit a first post-Newtonian expansion [21]. The techniques 
used in [20] and [21] are based on energy estimates derived from a first order nonlocal 
formulation of the Einstein-Euler equations in the harmonic gauge. A further improvement 
was achieved in [22] through the use of dispersive wave estimates in addition to the energy 
estimates. This improvement allowed the conditions on the initial data need to ensure the 
existence of a Newtonian limit to be relaxed as compared to [20] • It is also possible to use 
the dispersive wave estimates to extend the results of [22], and establish the existence of 
post-Newtonian expansions to the second order. 

The main aim on this paper is to describe, in a non-technical fashion, the rigorous results 
obtained in JSOJ [211 1221 [23] on the Newtonian limit and post-Newtonian expansions, and 
moreover, to relate these results to frame theory. These rigorous results, which originate 
from ideas introduced by Jiirgen, show that the question about the relationship between 
Einstein's and Newton's theories of gravity is not just philosophical, but in fact, leads to a 
deeper insight into the structure of the equations, and a better understanding of both the 
Newtonian limit and post-Newtonian approximations. We devote this article to Jiirgen's 
memory. 

2 The Frame Theory 

In this section, we briefly describe frame theory closely following Jiirgen's paper [6]. Frame 
theory is formulated using the following collections of fields on a spacetime (i.e. a 4- 
manifold) : 

tik a nowhere vanishing symmetric, 2-covariant tensor field, 

s tk a nowhere vanishing symmetric, 2-contravariant tensor field, 

Y\ k a symmetric linear connection, and 

T lk a symmetric, 2-contravariant tensor field. 

The fields t ik , s ik , Y\ k , and T ik are referred to as the the temporal metric, inverse temporal 
metric, gravitational field, and matter tensor, respectively. The laws of frame theory contain 
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these fields and two real-valued parameters A and A. To formulate the laws, it is convenient 
"raise" and "lower" tensorial indices using s lk and t ik , respectively. For example, we define 

V':=t ki V\ and := s ik u k . 

To avoid confusion it is necessary to indicate by dots whether the original object has upper 
or lower indices, as shown. 

The (local) laws of frame theory are: 

(i) For any timelike vector V 1 (i.e. ti k V l V k > 0), the quadratic form uji — > s lk 0JiU) k is 
positive definite on the subspace { Ui \ uOiV 1 = } of the covector space. 

(ii) The temporal and inverse temporal metrics are related to each other by 

t lk s kl = -XSl 

and satisfy 

t ik -i = 0, and s ik -i = 0, 
where the semicolon denotes the covariant derivative determined by the connection 

1 ik- 

(iii) The curvature tensor R? ikj of the connection Y l ik satisfies 

pi / _ pZ i 

while the Ricci tensor R ik = R? ' ik j and the matter tensor are related by 

Rik = 8n(T; k ' - f ik T j ]) - At ik . 

(iv) The matter tensor satisfies 

T k , k = 0. 

(v) In addition to general laws (i)-(iv), a matter model must be specified. Here, we consider 
only the case of a perfect fluid where 

T ik = (p + Ap)£P£/ fe + ps ik 

with t ik U i U k = 1, p > 0, and p + \p > 0. 
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The following facts can be inferred from the references [3... 9] in [6]. If A > 0, then, by 
appropriate rescaling, one can choose A = l,<7j& = —Uk, g lh = s lh , and in this case frame 
theory reduces to general relativity with the natural units G — c — 1, and a spacelike 
signature (+ + H — ). For the purpose of this paper, however, it is preferable not to fix the 
scaling and work with = —Xgik, s tk = g ik , and arbitrary A > 0. 

If A = 0, then frame theory reduces to Newton-Cartan theory. In this case, the temporal 
metric can be expressed in terms of a scalar absolute time t by 

tik t,it,ki 

and it is possible to choose (local) coordinates (x , x 4 = t) K = 1, 2, 3 such that 

t ik = diag(0, 0, 0, 1), s ih = diag(l, 1, 1, 0), 

and the nonzero components of the connection Y l kl are given by 

T? t =-g K and l£ = s KF e LFJ u J . 

Here, the spatial vector fields g(x,t) = (g K (x,t)) and u(x,i) = (u K (x,t)) play the part of 
the gravitational acceleration and the Coriolis angular velocity, respectively. While Newton- 
Cartan theory is more general than Newtonian gravity, it does reduces to Newton's theory 
if uj only depends on t. This is due to the existence of non-rotating coordinates for which 
uj = 0, and 

r kl = t tk t,is ij ^ tj , 

where $ is the Newtonian potential. 

The simplest condition that guarantees that Q depends only on time is due to Trautmann 
[3T] . and it reads 

R ik .lm = (2.1) 

or equivalently 

t,i[R J kl]m = . (2.2) 

This equation says that parallel transport of spacelike vectors V 1 , which are characterized 
by tjV* = 0, is path independent, a well known property of Newton's theory. If one restricts 
the solutions of the frame theory to spatially asymptotically flat ones (i.e. isolated systems), 
then for A = 0, the conditions fl2.ip - fl2.2l) are automatically satisfied, and Newtonian gravity 
is obtained. 

What is the physical meaning of A ? One interpretation is to consider all quantities 
having physical dimensions. Then A can be identified with 1/c 2 . To avoid confusions about 
the meaning of c — > oo, it is preferable to express the freedom of choosing units as certain 



5 



similarity transformations as discussed in [BJ. With this in mind, the Newtonian limit of a 
family 

{g ik (\,xi),T ik (\,xi) , A > 0} 

of solutions of Einstein's theory can be defined provided that g lk (\,x^), T lk (X,x^) and t^ = 
— \gik{\, x^) have, together with a number of derivatives, a limit for A = 0, and these limits 
satisfy the axioms (i)-(v) above for frame theory. Furthermore, if the family is differentiable 
in A, then one can obtain post-Newtonian expansions. 

3 The Newtonian limit 

In this section, we describe the approach taken in [2UJ ED [22] to mathematically justify the 
Newtonian limit for isolated systems. The articles [201 I2TI 122] cover the situation where the 
matter model is a perfect fluid, and so, this is what we restrict ourselves to in the following 
sections. However, it should be noted that the results of [20], [2TJ, [22] do not rely critically 
on the matter model, and can be extended to any other matter model that is well-posed, 
admits solutions with compact support, and has a matter propagation speed that remains 
bounded as c — > oo. Vlasov matter is an example of a model that satisfies these conditions. 

As shown in [20], the Einstein- Euler equations, which govern a gravitating perfect fluid, 
can be written in the following rescaled form 

G ij = 2e ±T ij and ViT ij = 0, (3.1) 

where 

T ij = (p + e 2 p)v i v j + pg ij and v\ = — -. 

e z 

In this formulation, the fluid four- velocity v\ the fluid density p, the fluid pressure p, the 
metric g^, and the coordinates (x l ) (i = 1, . . . ,4) are dimensionless. Spacetime M is taken 
to be diffeomorphic to IR 3 x [0, T) which is appropriate for modeling an isolated system. 
The coordinates (x 1 ) are global Cartesian coordinates on M = M. 3 x [0,T) with the (x 1 ) 
(I = 1,2,3) defining spatial coordinates on R 3 while t = x 4 is a Newtonian time coordinate 
on [0,T). 

As discussed in the introduction, the study of the equations (13.1 j) in the limit e \ is 
called the Newtonian limit. The overwhelming majority of investigations of the Newtonian 
limit, for example see [U El HU 15] and reference cited therein, have relied on formal ex- 
pansion in the parameter e. The main goal of the formal expansion is to identify suitable 
gravitational-matter variables and a gauge so that the gravitational-matter variables are well 
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defined in the limit e \ 0, and the Einstein-Euler equation (13.11) reduce to the Poisson-Euler 
equations 

dtp + d^/yw 1 ) = 0, (I, J =1,2, 3) (3.2) 

p(d t w J + w I drw J ) = -(p<9 J $ + d J p) , (d 1 = 5 IJ dj) (3.3) 

A$ = p, (A = 9 / 9 / ) (3.4) 

of Newtonian gravity in the limit e \ 0. However, as remarked in the introduction, this is far 
from what one actually wants, which is to be assured of the existence of a sufficiently large 
set of one-parameter families of solutions depending on e to the Einstein-Euler equations 
that exists on a common piece of spacetime of the form M = R 3 x [0, T), and converge, in a 
suitable sense, as e \ to a solution of the Poisson-Euler equations on M. The sufficiently 
large set of one-parameter families of solutions should be large enough so that there are not 
any unreasonable restrictions on the choice of the initial data for the matter, which in the 
perfect fluid case is the fluid density and spatial three-velocity. If this can be achieved, then 
this is what is meant by a rigorous Newtonian limit. 

The rigorous analysis of the Newtonian limit starts by viewing the Einstein-Euler equa- 
tions as an initial value problem for a gauge hyperbolic system. The main difficulty in 
obtaining a rigorous Newtonian limit is that the limit e \ is a singular limit for the 
Einstein-Euler equations. Roughly speaking, this means that terms of the form 1/e will 
appear in the evolution equations, and these terms become unbounded as e \ 0. This leads 
to the following mathematical problems that must be sequentially solved. 

(i) The first problem is to find a one parameter family of initial data 

{fl'y|t=o,^fl'y|t=o J Pe|t=o,^|t=o} < e < e (3.5) 
on the spacelike hypersurface 

Z = {(x I ,0)\(x I )eR 3 } 

that solves the constraint equations, and converge, in an appropriate sense, as e \ 
to Newtonian initial data. 

(ii) Once problem (i) is solved and an appropriate gauge is chosen, local existence theorems 
for hyperbolic partial differential equations show that the one parameter family of 
initial data (13.51) generates a one-parameter family of solutions 

{gUx, t),p £ (x, t),vl(x, t)} < e < e (3.6) 
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on a spacetime region of the form 

(x,t) e M e = M 3 x [0,T £ ). 

The second problem is to show that there exists aT, > such that T t > T* for all 
< e < eo, because otherwise the one parameter family of solutions (13.61) would not 
exist on a common piece of spacetime, and consequently, it would be impossible to 
even discuss the limit e \ 0. 

(iii) A solution to problem (ii) guarantees the existence of a common piece of spacetime 
M* = [0, T*) x IR 3 on which the one-parameter family of solutions (13. 6p to the Einstein- 
Euler equations exist. So once problem (ii) is solved, the final problem is to demonstrate 
the convergence on M* as e \ of the Einstein-Euler solutions (13.61) to solutions of 
the Poisson-Euler equations. This convergence should be uniform in e, and thus, part 
of the problem is to show the existence of an order e error estimate that measure the 
difference between the fully relativistic and Newtonian solutions. 

A positive resolution to these three questions is the minimum requirement for a rigorous 
Newtonian limit. More information about the limit e \ would certainly be desirable, and, 
in some cases, essential for answering questions about the behaviour of the Newtonian limit. 
For instance, rather than a lower bound T*, it would be highly desirable to determine the 
maximal interval T(e) > T* on which the family of solutions (13.61) exist, and converge to 
solutions of the Poisson-Euler equations. However, this is an extremely difficult problem that 
is at least as difficult as determining the long time behavior of the Poisson-Euler equations, 
which is still an unresolved problem. 

There have been various methods proposed to solve the three problems listed above. 
Each method is dependent on the gauge and the gravitational-matter variables that are 
chosen. As stated above, we will outline the method used in [2D]. For different approaches, 
see J2DJ [TDJ . The analysis in [2D] starts by replacing the metric g^ and the fluid velocity v l 
with new variables that are compatible with the limit e \ 0. The new gravitational variable 
is a density defined via the formula 

9 ij = , 6 Q lj (3.7) 
y/-det{Q) 

where 

™ f$ IJ 0\ 2 (Au IJ \ „ 3 / u /4 \ , 4 fO \ 

From this, it not difficult to see that the density u JJ is equivalent to the metric gij for e > 0, 
and is well defined at e = 0. For the fluid, a new velocity variable w l is defined by 

v 4 — 1 

v 1 = w 1 and w 4 = . (3.9) 
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The choice of the variables u ,J for the gravitational field was inspired by the gravitational 
variables discovered by Jiirgen, see Figured] at the end of the article. To see the relationship, 
we let u l J denote Jiirgen's variables. In terms of these variables, the metric is given by the 
formula 

where 



e ~e 



The formulas (13. 7p . (13.81) . (I3.10p . and (13. lip show that u tj are related to the u l J through a 
simple rescaling of the components by an appropriate power of e. 
For technical reasons, an isentropic equation of state 

p = Kp (n+1)/n , (3.12) 

where K e R>o and n G N is assumed. This allows for the use of a technique of Makino [T7j 
to regularize the fluid equations through the use of the following fluid density variable 

^ ( 4A> i( » + i)r a2 "- (3 - 13) 

This choice makes it possible to formulate the Euler equations as a hyperbolic system that 
is regular across the fluid-vacuum interface. In this way, it is possible to construct solutions 
to the Einstein-Euler equations that represent compact gravitating fluid bodies (i.e. stars) 
both in the Newtonian and relativistic setting [T71 [23] . In the Newtonian setting, this is 
straightforward to see. Using (13.121) and (13.131) . the Poisson- Euler equations (I3.2l) - (l3.4p 
become 

o 

d t a = —w I dia — —diw 1 , (3-14) 

d tW J = d J a - w^rW 1 - <9 J $, (3.15) 

2n 

A$ = p, (3.16) 

which is readily seen to be regular even across regions where a vanishes. 

It is important to note that the Makino technique only produces a particular subclass 
of solutions to the Euler equations representing compact fluid bodies. In particular, the 
type of fluid solutions obtained by this method have freely falling boundaries [25], and 
hence, do not include static stars of finite radius. On the one hand, this makes the Makino 
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technique far from ideal, but on the other, it is the only method presently available for 
generating relativistic solutions that represent gravitating, compact fluid bodies. However, 
in trying to understand the Newtonian limit and post-Newtonian approximations, these 
solutions are almost certainly general enough to obtain a comprehensive understanding of the 
mathematical issues involved in the Newtonian limit and post-Newtonian approximations. 
This is because the mathematical difficulties in studying the post-Newtonian expansions 
do not rely heavily on the structure of the matter model. As stated above, the matter 
model need only be well-posed, admit solutions with compact support, and have a matter 
propagation speed that remains bounded in the limit e \ 0. 

The approach taken in [201 (2TJ [22] relies on a particular nonlocal first order formulation 
for the Einstein-Euler system. This formulation starts by the introduction of the following 
vector 

V = (edtu^ d/uVu^', a, vf), 

which contains both the gravitational and matter variables. Assuming a harmonic gauge, a 
straightforward calculation shows that V satisfies an equation of the form 

b°(eV)d t V = -c'djV + &'(e, V)djV + f(e, V) + -g(p) (3.17) 
e e 

where b°(eV) and 6 / (e, V) are symmetric matrices that depend analytically on eV and (e, V), 
respectively, / is a map that depends analytically on (e, V), 

(A 

5 IJ 

0, 



0/ 

and 

g= 0,0, 0,0) T . 

This formulation has the virtue that the are only two terms which contain 1/e, and these 
terms are of a particularly simple structure. However, from an analysis point of view, 
the presence of the term -g(p) is problematic. Fortunately, this term can be removed by 
subtracting off the gravitational contribution from the Newtonian potential in the following 
fashion 

W = (ed t u ij , dju ij - <5|5id/$, eu ij , a, w*) (3.18) 



(5 1J 
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where $ solves 

A$ = p. 

This effectively moves the problem term -g(p) into the c l diW term. In terms of the new 
variable W, the Einstein-Euler equations become 

b°(eW)d t W = - JdjW + 6 J (e, W)djW e + F(e, W), (3.19) 
e 

where F is a nonlocal map that depends analytically on (e, W). The importance of this non- 
local symmetric hyperbolic formulation of the Einstein-Euler equations is that it is possible, 
given suitable initial data, to use energy estimates to establish the existence of solutions on 
a time interval that is independent of e. In fact, equation (I3.19P is a particular example 
of a class of singular hyperbolic equations for which the limit e \ has been well studied 
[TTI [T2"| I2"8"| [291 150] . However, it is important to note the results of these articles do not 
directly apply equation (13.191) . The reason for this is that the initial data for the system 
(I3.19P must include a 1/r piece for the metric, and hence, cannot lie in the Sobolev space 
H k on which the results of pH QH [281 ISS [30] are based. 

As discussed above, the first problem is to generate suitable initial data for the evolution 
equation ( I3.19j) . The method used in [20], [22] for generating this data is based on a method 
introduced by Lottermoser [16] . The result is a positive solution to problem (i) listed above. 
We will not discuss the details here and only outline the main result. 

First, let {5, s, R, k} be a fixed set of numbers satisfying — 1 < 5 < — 1/2, s G Z> 2 , R > 0, 
k G Z> 3+s . Then it is shown in [22J, that for any choice of the free initial data 

(/ J , ii J G Hi,, ao, wi) e H k 5 +1 x x x (3.20) 

where the density has compact support supp «o C Br, there exists a one parameter family 
of determined data 

to* O G tf, fc+1 x H k 5 +l x e G (-e , e ) (3.21) 

that depends analytically on e. The gravitational variable u l J and its first time derivative 
are determined on the hypersurface S in terms of the free and determined initial data as 
follows: 

(u«)=( S " J), (3.22) 



and 



IJ fi KI S 
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Together, <KT3f . <Km . and 

( Po ,wlKe) ( 3 - 24 ) 
solve the gravitational constraint equations G I0 —2e 4 T 10 = 0, the harmonic condition g^T^ = 
0, and the fluid velocity normalization v % v^ = — 1/e 2 on S for — eo < e < eo- Moreover, </> e , 
rof , and Wq € satisfy the equations 

A0 o = p o + 9 7 V J , An^ = -<9 L ^ J , and w£ = , 

respectively, at e = 0. We also note that the free data ( 13.201) can be chosen so that any one 
of the 3 /J , 34 J , a , and depends analytically on e. Indeed, as will be discussed below, the 
rigorous post-Newtonian expansions require that the gravitational free initial data $ IJ , $ I 4 J 
is e-dependent. 

With the question of generating suitable initial data resolved, the next step is to resolve 
problem (ii) which concerns finding a lower bound on the time of existence for solutions 
W € (x,t) to 03.191) generated from the e-dependent initial data fl3.22p - fl3.24p . This requires 
decomposing W e into the sum of a time independent 1/r term and a time varying 1/r 2 term: 

W e (x,t) = W- lte (x) + W- 2 , e (x,t), 

where 

WC 1 (a;) = (0,0,e< J '| t= o,0,0) and W- 3>e (x,t) = W e (x,t) -W-^x). 
The initial data given by fl3.22p - fl3.24p then implies that 

||W_ 1)e |U + i < 1 and ||W_ 2 , e (0)|U <1 VeG[0,e ). 

5 o—l 

This bound on the initial data and the special form of equation (13.191) are the key to deriving 
e-independent energy estimates. These energy estimates provide a uniform bound on the 
solution for some time T* which may, a priori, be larger than the time of existence. But, 
by continuation principle for symmetric hyperbolic equations, solutions can be continued as 
long as they remain bounded in the W 1,oc norm, and so, the e-independent energy estimates 
together with the Sobolev inequality guarantee the existence of a classical solution 

W e (x,t) G C 1 (IR 3 x [0, T*)) 0<e<e 

to (13.191) that satisfies the uniform bound 

H j (t)he + \\W- 2 M\h» + e\\dtW e {t)\\ H *-i < 1 Ve e (0,e ). 

This not only resolves problem (ii) stated above, but also provides uniform estimates for the 
one parameter family of solutions to the Einstein- Euler equations determined by W e (x,t). 
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For initial data that is general as fl3.22p - fl3.24p . energy estimates alone are not sufficient 
to resolve problem (iii), and establish convergence to a Newtonian solution. To obtain 
convergence, the energy estimates must be combined with dispersive estimates for the wave 
equation [22]. Together, these estimates can be used to establish the existence of a solution 
{p(x,t),w I (x,t),^(x,t)} to the Poisson-Euler equations 

d t p + d I (pw I )=0, 
p(d t w J + w 7 9/w J ) = -(p<9 J i> + d J p) , 
A$ = p 

with initial conditions 

p\ t=0 = {AKn{n + l)) -1 ^™, u> 7 | i= o = ™o 
and a solution u e (x,t) of the wave equation 

e 2 d^ - Afi« = -SiSip + e 2 6idid^, (3.25) 

with initial conditions 

dtK 3 \ t= , = \{^Wa - 2d l2 IJ s { ;5^ + s i AA- 1 & IJ tL J ) + 5m$\ t=0 , 

such that the following uniform estimates hold 

H j (t) - u«|| z « + \\dMt) - 9/U«(t)|| H *- a + \\ed t ui j (t) - ed#i{t)\\ Hk -i < e (3.26) 

and 

\\Pe{t) -p(t) \\ H k-2 + \\wiit) -w T (t) \\ Hk - 2 + K(t)|| H *- a < e, (3.27) 

for all (t,e) G [0, T*) x (0,e ). It is this last estimate (13.261) that resolves problem (iii) 
and established rigorously the convergence as e \ of the fully relativistic solution of the 
Einstein-Euler equations determined by W e (x, t) to a solution of the Poisson-Euler equations. 
Clearly, it also provides an order e error estimate that measures the difference between the 
relativistic and Newtonian solutions. We also note that 

||^'IU« + \\d&*(t)\\H>-* <1 (3.28) 

for t e [0, T ). Combined, the formula ()3.7|) . and the estimates (I3.26P and (13.281) show that 

H 3 (t) - SV&'SuWv + \\d ig v(t)\\ Hk - 2 < e (3.29) 
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and 

\\e 2 gl 3 {t) + 6t6f\\ L e + W^djgt^Wn^ < e (3.30) 

for all t G [0,T ). It also follows from the estimates (I3.26p - fl3.28p . and the fact the support 
of p e stay uniformly bounded, that the stress energy tensor T % j satisfies 

\\T?(t) - ^{P{tW{tW{t)+Pm IJ )\\H^ < e (3.31) 

for all t G [T , e). 

The estimate (13.261) is not strong enough to guarantee that the gravitational variables 
u u converge to a particular limit. It only says the gravitational part of the fully relativistic 
solution is approximated to order e by a solution of the singular wave equation (13.251) . 
However, after an arbitrarily small time, the dispersive effects of the wave equation are 
strong enough to allow for convergence, at least for the spatial and temporal derivatives. 
This can be seen by the following estimate: given any t G (0,T,), the solution u* J satisfies 

- 5l5idMt)\\ wk -2^ + e\\d£ ij (t) - 5i5idMt)\\ wk -2^ < t0 e (3.32) 

for all t G (to, T*). Moreover, restricting to any ball B\ C K 3 of radius A, the above estimate 
improves to 

\\d^{t) - 5l5ldMt)\\w*-^ <t ,A e 3/2 , (3.33) 

and 

\\d t U ij (t) - 5\5{d t ^{t)\\ w u-,,o a{BA) < t0)A (3.34) 

for all t G (t ,T*). We also note that in the estimates fl3~26D - fl3T27D . and fl3~32|) - (l3T3D . each 
time a time derivative is taken, the estimate gets worse by a power of e. For example, it 
follows from (13.341) that 

\\d^(t) - 8i8id 2 Mt)\\ w u-^ {BA) < t0 ,A \ (3.35) 

for all t G (to, T*) 

Together, the estimates ( 13.261) . and (I3.32p -( l3~34l) establish that diu l J and d t u l J converge 
to 5l5{di& and 5l8{d t &, respectively, as e \ 0, at least after waiting a short time. They also 
show that the Christoffel symbols satisfy 

||r* - %fi}5 Ik dMt)\\ w >s.- < t0 e (3.36) 

and 

HftT* - 5t5p Ik djdMt)\\w^(B A) <t ,A (3.37) 
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for all t G (to,T m ). In particular, this implies that on the spacetime region M\ tto = x 
(to,T*), both the Christoffel symbols and the curvature components converge as e \ 0, and 
hence this limit satisfies the assumptions on Jiirgen's frame theory on the bounded spacetime 
region M A)V 

It is possible to prove convergence of the relativistic solution to a Newtonian solution 
using energy estimates alone if the initial data is suitably restricted. Indeed, as shown in 
|2Uj . if the free initial data is chosen so that 

3 " = 0(e) and i\ 3 = 0(e) as e \ 0, 

then W e (t, x) satisfies 

\\d t W e \ t= 4 Hk -,<l ee(0,e ). (3.38) 

This uniform estimate combined with energy estimates can be used to improve (I3.26P and 
(153711 to 

\\p e (t) - p{t)\\ Hk -i + K(t) -w I {t)\\ H i.-i + K 4 (i)||^-i < e (3.39) 

and 

\\v?i{t) - 5{5{m\y + ||0,u*(*) - 5\5{dMt)\\H^ + ||e^'|| H *-i < e. (3.40) 

In particular, this shows that the condition (I3.38P on the initial data implies that both the 
gravitational u l J and the matter variables {p e ,w l e } converge as e \ 0. We also note that 
with this restricted data, the estimates (13.361) and (I3.37P improve to 

\\r^ - stsp^dAm^ < e 

for all t e [0, T*) and 

||^r* - 5t5*5 Ik d I dMt)\\w*-^ < t0 ,A ^ 

for all t G (to, T*), respectively. This gives convergence on the space time region M to = 
M 3 x(t ,T,). 

To get convergence, in the sense of frame theory, on the whole spacetime slab M* = 
M 3 x [0,T*), requires that initial data is chosen so that both the first and second time 
derivatives of W e at t — are bounded as e \ 0. This type of additional restriction of the 
time derivatives is also what is needed to obtain rigorous post-Newtonian expansions, which 
will be discussed in the next section. 
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4 Post-Newtonian expansions 



That the boundedness condition (13.381) on the first time derivative of W e at t = implied 
the existence of a th order expansions in e is a manifestation of Kreiss 's bounded derivative 
principle [2, [13]. This principle states that using energy estimates alone, the problem of 
generating £ th order expansions in e for solutions to an equation of the form CI3. 19f) can be 
reduced to the problem of finding initial data that satisfies 

II Wt=o|U- f < 1 (4.1) 

o — 1 

for p = 1, 2, ...,£+ 1. This process of choosing initial data to satisfy (14.11) is called initial- 
ization. This shows that the existence problem for the post-Newtonian expansions can be 
replaced by the problem of finding initial data that can be properly initialized. We note 
that the conditions (14. Xp effectively amounts to additional elliptic equations that the initial 
data must satisfy in addition to the usual constraint equations. For an example of these 
additional elliptic equations on the initial data, see equation (6.4) in [2~T] . 

One of the results of [21] is to show that it is possible to construct initial data that 
satisfies 

\\%W e \ t=0 \\ Hk - P < 1 p = 0,l,2,3 (4.2) 

5 — 1 

By Kreiss's bounded derivative principle, this is enough to prove the existence of a first 
post-Newtonian expansion for the Einstein-Euler equations. More specifically, it is possible 
to show that the solution 

K J '(x, t), a e (x, t), w\{x, t)} (x, t) e M* 

admit convergent expansions (uniform for < e < eo) of the form 

o 2 00 

q=l q=3 

2 °° 

q=\ q=3 
2 °° 

e v d v t v?i = e v 8\8id v t <& + ^ e q+u d^u ij + ^ e q+v d^uf 

q=l q=3 



(4.3) 

u = 0,l, (4.4) 
^ = 1,2, (4.5) 
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and 

2 oo 

d v t a t = dfa + e<ld t® + Yl e<?(9 ^ v = 0,l, (4.6) 

g=l q=3 
2 oo 

g=l <?=3 

where the first expansion is convergent in C°([0, T*);L|), and the rest are convergent in 

C°([0, T*); H k ~^). Here, the expansions coefficient (g = 1,2) {u u (t), a(t), satisfies a 

linear (nonlocal) symmetric hyperbolic system that only depends on {a(t),w I (t), <&(£)} if 5 = 

1, and {ait),^ 1 ^),^),^^),^),^^)} if g = 2, while for g e Z> 3 , {%'(*). w l £ (t)} 
satisfies a linear (non-local) symmetric hyperbolic system that only depends on e, {a(t), w^t), $ 

a(f), &*(*)} for p= 1,2, and {u/ j (t), a e (t), &»(*)} for p = 3, 4, . . . , q - 1. 

Since, the {u lJ '(£), a(t), w l (t)} are e independent, the expansions fl4.3p - fl4.7l) do, in fact, 
coincide with the first post-Newtonian expansion. This can be seen by introducing the 
following e-independent quantities: 

h ij = (4^'-2 % i fc V) 9 = 1,2, 

where (77^) = diag(l, 1, 1, — 1). From these expression and the expansions (14.31) - (14. 7p . it 
follows that (see (13. 7p - (13.81) ) the metric gij is given by 

9u = -1 - 2$ - e/> 44 - 6 2 (3(!>) 2 + + 0(e 3 ), 

and 

g rj = Su - 2e 2 5 u l> - e 3 h IJ - e 4 ((£)%./ + h IJ ) + 0(e 5 ). 

It is important to note that higher order expansions in e can be generated for the metric 
Qij using fl4.3l) - fl4.7p . However, these higher order terms will, in general, depend on e in a 
non-analytic fashion. Consequently, they represent some sort of generalized post-Newtonian 
expansion, and it is not clear without further analysis how these expansions are related to 
the standard post-Newtonian ones. 
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To make progress beyond the first post-Newtonian expansion, the simplest thing to do 
would be to try and prove the existence of initial data that satisfies (14.1 1) for £ > 3. In 
the present setup using the harmonic gauge, the condition (14. ip for £ = 4 leads to elliptic 
equations that are not solvable within the required function spaces. On the other hand, it is 
expected that with a suitable gauge choice, it should be possible to generate post-Newtonian 
expansions to at least the 2.5 post-Newtonian order after which there are indications that 
the post-Newtonian expansions will break down. For a lucid discussion of this phenomenon 
see [21]. 

As remarked in [24] , the choice of harmonic gauge may be the reason for not being able 
to reach the 2.5 post-Newtonian order, and there may well exist other gauges for which the 
conditions (14.1 j) can be satisfied for £ > 3. However, even if this is the case, it would still need 
to be verified that these other gauges are compatible with the singular hyperbolic energy 
estimates that are guaranteed to arise. 

Kreiss's bounded derivative principle which reduces the existence problem of £ th order 
expansions in e to that of finding initial data satisfying (14. ip relies on energy estimates. In our 
situation, we also have dispersive wave estimates coming from the fact that the gravitational 
variables satisfy a wave equation. These dispersive estimates can be used to show that 
the matter variables are more regular in e as compared to the gravitational ones. This means 
that it is possible to go beyond the first post-Newtonian expansion while only satisfying (14. ip 
for £ = 3. This type of effect was already established in [22] at the Newtonian level. This 
is clear from the estimates (I3.26I) - (I3.27I) which establish the existence of a Newtonian limit 
without need for the condition (13.381) . In fact, it is possible to show using the results of 
both the paper [2TJ [22] that there exists post-Newtonian expansions to the second order. 
However, it still remains to be seen what is the optimal choice of gauge that allows one to 
satisfy (14. ip to as high order as possible while still being able to use energy and dispersive 
estimates to generate post-Newtonian expansions to the highest possible order. 

It is well known that Newtonian gravity generalizes to the cosmological setting [27] . Using 
similar techniques as outlined above, it is also possible to rigorously establish the existence 
of cosmological post-Newtonian expansions [23] . In contrast to the asymptotically flat case, 
it is possible to construct initial data that satisfies (14. ip for arbitrarily large £. This implies 
that there is no obstruction to generating rigorous post-Newtonian expansions to arbitrary 
order which is somewhat surprising in light of difficulties that are encountered at both the 
rigorous and formal level in trying to develop post-Newtonian expansions on asymptotically 
flat spacetimes beyond the order 2.5. On asymptotically flat spacetimes, the problems that 
occur in the higher order post-Newtonian expansions are often attributed to the reaction 
of gravitational radiation with itself and matter. The analysis contained in the paper [23] 
shows that this is not the complete story as these effects are also present in the cosmological 
setting but do not cause similar difficulties. 
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\ ■ * 

Figure 1: Jiirgen's notes where he introduces the gravitational variables fl3.10l) - fl3.11l) . 
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